Hox genes and patterning the vertebrate body.

The diversity of vertebrate body plans is dizzying, yet stunning for the many things they have in common. Vertebrates have inhabited virtually every part of the earth from its coldest to warmest climates. They locomote by swimming, flying, walking, slithering, or climbing, or combinations of these behaviors. And they exist in many different sizes, from the smallest of frogs, fish and lizards to giraffes, elephants, and blue whales. Despite these differences, vertebrates follow a remarkably similar blueprint for the establishment of their body plan. Within the relatively small amount of time required to complete gastrulation, the process through which the three germ layers, ectoderm, mesoderm, and endoderm are created, the embryo also generates its body axis and is simultaneously patterned. For the length of this axis, the genes that distinguish the neck from the rib cage or the trunk from the sacrum are the Hox genes. In vertebrates, there was evolutionary pressure to maintain this set of genes in the organism. Over the past decades, much has been learned regarding the regulatory mechanisms that ensure the appropriate expression of these genes along the main body axes. Genetic functions continue to be explored though much has been learned. Much less has been discerned on the identity of co-factors used by Hox proteins for the specificity of transcriptional regulation or what downstream targets and pathways are critical for patterning events, though there are notable exceptions. Current work in the field is demonstrating that Hox genes continue to function in many organs long after directing early patterning events. It is hopeful continued research will shed light on remaining questions regarding mechanisms used by this important and conserved set of transcriptional regulators.

A novel approach to exploring the dark genome and its application to mapping of the vertebrate virus fossil record.

BACKGROUND: Genomic regions that remain poorly understood, often referred to as the dark genome, contain a variety of functionally relevant and biologically informative features. These include endogenous viral elements (EVEs)-virus-derived sequences that can dramatically impact host biology and serve as a virus fossil record. In this study, we introduce a database-integrated genome screening (DIGS) approach to investigate the dark genome in silico, focusing on EVEs found within vertebrate genomes. RESULTS: Using DIGS on 874 vertebrate genomes, we uncover approximately 1.1 million EVE sequences, with over 99% originating from endogenous retroviruses or transposable elements that contain EVE DNA. We show that the remaining 6038 sequences represent over a thousand distinct horizontal gene transfer events across 10 virus families, including some that have not previously been reported as EVEs. We explore the genomic and phylogenetic characteristics of non-retroviral EVEs and determine their rates of acquisition during vertebrate evolution. Our study uncovers novel virus diversity, broadens knowledge of virus distribution among vertebrate hosts, and provides new insights into the ecology and evolution of vertebrate viruses. CONCLUSIONS: We comprehensively catalog and analyze EVEs within 874 vertebrate genomes, shedding light on the distribution, diversity, and long-term evolution of viruses and reveal their extensive impact on vertebrate genome evolution. Our results demonstrate the power of linking a relational database management system to a similarity search-based screening pipeline for in silico exploration of the dark genome.

Self-supervised learning on millions of primary RNA sequences from 72 vertebrates improves sequence-based RNA splicing prediction.

Language models pretrained by self-supervised learning (SSL) have been widely utilized to study protein sequences, while few models were developed for genomic sequences and were limited to single species. Due to the lack of genomes from different species, these models cannot effectively leverage evolutionary information. In this study, we have developed SpliceBERT, a language model pretrained on primary ribonucleic acids (RNA) sequences from 72 vertebrates by masked language modeling, and applied it to sequence-based modeling of RNA splicing. Pretraining SpliceBERT on diverse species enables effective identification of evolutionarily conserved elements. Meanwhile, the learned hidden states and attention weights can characterize the biological properties of splice sites. As a result, SpliceBERT was shown effective on several downstream tasks: zero-shot prediction of variant effects on splicing, prediction of branchpoints in humans, and cross-species prediction of splice sites. Our study highlighted the importance of pretraining genomic language models on a diverse range of species and suggested that SSL is a promising approach to enhance our understanding of the regulatory logic underlying genomic sequences.

Paired fins in vertebrate evolution and ontogeny.

The origin of paired appendages became one of the most important adaptations of vertebrates, allowing them to lead active lifestyles and explore a wide range of ecological niches. The basic form of paired appendages in evolution is the fins of fishes. The problem of paired appendages has attracted the attention of researchers for more than 150 years. During this time, a number of theories have been proposed, mainly based on morphological data, two of which, the Balfour-Thacher-Mivart lateral fold theory and Gegenbaur's gill arch theory, have not lost their relevance. So far, however, none of the proposed ideas has been supported by decisive evidence. The study of the evolutionary history of the appearance and development of paired appendages lies at the intersection of several disciplines and involves the synthesis of paleontological, morphological, embryological, and genetic data. In this review, we attempt to summarize and discuss the results accumulated in these fields and to analyze the theories put forward regarding the prerequisites and mechanisms that gave rise to paired fins and limbs in vertebrates.

Shark genome size evolution and its relationship with cellular, life-history, ecological, and diversity traits.

Among vertebrates, sharks exhibit both large and heterogeneous genome sizes ranging from 2.86 to 17.05 pg. Aiming for a better understanding of the patterns and causalities of shark genome size evolution, we applied phylogenetic comparative methods to published genome-size estimates for 71 species representing the main phylogenetic lineages, life-histories and ecological traits. The sixfold range of genome size variation was strongly traceable throughout the phylogeny, with a major expansion preceding shark diversification during the late Paleozoic and an ancestral state (6.33 pg) close to the present-day average (6.72 pg). Subsequent deviations from this average occurred at higher rates in squalomorph than in galeomorph sharks and were unconnected to evolutionary changes in the karyotype architecture, which were dominated by descending disploidy events. Genome size was positively correlated with cell and nucleus sizes and negatively with metabolic rate. The metabolic constraints on increasing genome size also manifested at higher phenotypic scales, with large genomes associated with slow lifestyles and purely marine waters. Moreover, large genome sizes were also linked to non-placental reproductive modes, which may entail metabolically less demanding embryological developments. Contrary to ray-finned fishes, large genome size was associated neither with the taxonomic diversity of affected clades nor with low genetic diversity.

In Silico Analyses of Vertebrate G-Protein-Coupled Receptor Fusions United With or Without an Additional Transmembrane Sequence Indicate Classification into Three Groups of Linkers.

Natural G-protein-coupled receptors (GPCRs) rarely have an additional transmembrane (TM) helix, such as an artificial TM-linker that can unite two class A GPCRs in tandem as a single-polypeptide chain (sc). Here, we report that three groups of TM-linkers exist in the intervening regions of natural GPCR fusions from vertebrates: (1) the original consensus (i.e., consensus 1) and consensus 2~4 (related to GPCR itself or its receptor-interacting proteins); (2) the consensus but GPCR-unrelated ones, 1~7; and (3) the inability to apply 1/2 that show no similarity to any other proteins. In silico analyses indicated that all natural GPCR fusions from Amphibia lack a TM-linker, and reptiles have no GPCR fusions; moreover, in either the GPCR-GPCR fusion or fusion protein of (GPCR monomer) and non-GPCR proteins from vertebrates, excluding tetrapods, i.e., so-called fishes, TM-linkers differ from previously reported mammalian and are avian sequences and are classified as Groups 2 and 3. Thus, previously reported TM-linkers were arranged: Consensus 1 is [T(I/A/P)(A/S)-(L/N)(I/W/L)(I/A/V)GL(L/G)(A/T)(S/L/G)(I/L)] first identified in invertebrate sea anemone Exaiptasia diaphana (LOC110241027) and (330-SPSFLCI-L-SLL-340) identified in a tropical bird Opisthocomus hoazin protein LOC104327099 (XP_009930279.1); GPCR-related consensus 2~4 are, respectively, (371-prlilyavfc fgtatg-386) in the desert woodrat Neotoma lepida A6R68_19462 (OBS78147.1), (363-lsipfcll yiaallgnfi llfvi-385) in Gavia stellate (red-throated loon) LOC104264164 (XP_009819412.1), and (479-ti vvvymivcvi glvgnflvmy viir-504) in a snailfish GPCR (TNN80062.1); In Mammals Neotoma lepida, Aves Erythrura gouldiae, and fishes protein (respectively, OBS83645.1, RLW13346.1 and KPP79779.1), the TM-linkers are Group 2. Here, we categorized, for the first time, natural TM-linkers as rare evolutionary events among all vertebrates.

Evolutionary origins of bitter taste receptors in jawed vertebrates.

Taste is a sense that detects information about nutrients and toxins in foods. Of the five basic taste qualities, bitterness is associated with the detection of potentially harmful substances like plant alkaloids. In bony vertebrates, type 2 taste receptors (T2Rs), which are G-protein-coupled receptors (GPCRs), act as bitter taste receptors1,2. In vertebrates, six GPCR gene families are described as chemosensory receptor genes, encoding taste receptor families (T1Rs and T2Rs) and olfactory receptor families (ORs, V1Rs, V2Rs, and TAARs). These families of receptors have been found in all major jawed vertebrate lineages, except for the T2Rs, which are confined to bony vertebrates3. Therefore, T2Rs are believed to have emerged later than the other chemosensory receptor genes in the bony vertebrate lineage. So far, only the genomes of two cartilaginous fish species have been mined for TAS2R genes, which encode T2Rs4. Here, we identified novel T2Rs in elasmobranchs, namely selachimorphs (sharks) and batoids (rays, skates, and their close relatives) by an exhaustive search covering diverse cartilaginous fishes. Using functional and mRNA expression analyses, we demonstrate that their T2Rs are expressed in the oral taste buds and contribute to the detection of bitter compounds. This finding indicates the early origin of T2Rs in the common ancestor of jawed vertebrates.

Deep mining of the Sequence Read Archive reveals major genetic innovations in coronaviruses and other nidoviruses of aquatic vertebrates.

Virus discovery by genomics and metagenomics empowered studies of viromes, facilitated characterization of pathogen epidemiology, and redefined our understanding of the natural genetic diversity of viruses with profound functional and structural implications. Here we employed a data-driven virus discovery approach that directly queries unprocessed sequencing data in a highly parallelized way and involves a targeted viral genome assembly strategy in a wide range of sequence similarity. By screening more than 269,000 datasets of numerous authors from the Sequence Read Archive and using two metrics that quantitatively assess assembly quality, we discovered 40 nidoviruses from six virus families whose members infect vertebrate hosts. They form 13 and 32 putative viral subfamilies and genera, respectively, and include 11 coronaviruses with bisegmented genomes from fishes and amphibians, a giant 36.1 kilobase coronavirus genome with a duplicated spike glycoprotein (S) gene, 11 tobaniviruses and 17 additional corona-, arteri-, cremega-, nanhypo- and nangoshaviruses. Genome segmentation emerged in a single evolutionary event in the monophyletic lineage encompassing the subfamily Pitovirinae. We recovered the bisegmented genome sequences of two coronaviruses from RNA samples of 69 infected fishes and validated the presence of poly(A) tails at both segments using 3'RACE PCR and subsequent Sanger sequencing. We report a genetic linkage between accessory and structural proteins whose phylogenetic relationships and evolutionary distances are incongruent with the phylogeny of replicase proteins. We rationalize these observations in a model of inter-family S recombination involving at least five ancestral corona- and tobaniviruses of aquatic hosts. In support of this model, we describe an individual fish co-infected with members from the families Coronaviridae and Tobaniviridae. Our results expand the scale of the known extraordinary evolutionary plasticity in nidoviral genome architecture and call for revisiting fundamentals of genome expression, virus particle biology, host range and ecology of vertebrate nidoviruses.

Genomic Resources and Annotations for a Colonial Ascidian, the Light-Bulb Sea Squirt Clavelina lepadiformis.

Ascidian embryos have been studied since the birth of experimental embryology at the end of the 19th century. They represent textbook examples of mosaic development characterized by a fast development with very few cells and invariant cleavage patterns and lineages. Ascidians belong to tunicates, the vertebrate sister group, and their study is essential to shed light on the emergence of vertebrates. Importantly, deciphering developmental gene regulatory networks has been carried out mostly in two of the three ascidian orders, Phlebobranchia and Stolidobranchia. To infer ancestral developmental programs in ascidians, it is thus essential to carry out molecular embryology in the third ascidian order, the Aplousobranchia. Here, we present genomic resources for the colonial aplousobranch Clavelina lepadiformis: a transcriptome produced from various embryonic stages, and an annotated genome. The assembly consists of 184 contigs making a total of 233.6 Mb with a N50 of 8.5 Mb and a L50 of 11. The 32,318 predicted genes capture 96.3% of BUSCO orthologs. We further show that these resources are suitable to study developmental gene expression and regulation in a comparative framework within ascidians. Additionally, they will prove valuable for evolutionary and ecological studies.

Conserved Noncoding Elements Evolve Around the Same Genes Throughout Metazoan Evolution.

Conserved noncoding elements (CNEs) are DNA sequences located outside of protein-coding genes that can remain under purifying selection for up to hundreds of millions of years. Studies in vertebrate genomes have revealed that most CNEs carry out regulatory functions. Notably, many of them are enhancers that control the expression of homeodomain transcription factors and other genes that play crucial roles in embryonic development. To further our knowledge of CNEs in other parts of the animal tree, we conducted a large-scale characterization of CNEs in more than 50 genomes from three of the main branches of the metazoan tree: Cnidaria, Mollusca, and Arthropoda. We identified hundreds of thousands of CNEs and reconstructed the temporal dynamics of their appearance in each lineage, as well as determining their spatial distribution across genomes. We show that CNEs evolve repeatedly around the same genes across the Metazoa, including around homeodomain genes and other transcription factors; they also evolve repeatedly around genes involved in neural development. We also show that transposons are a major source of CNEs, confirming previous observations from vertebrates and suggesting that they have played a major role in wiring developmental gene regulatory mechanisms since the dawn of animal evolution.

Comparative Analysis and Phylogenetic Study of Dawkinsia filamentosa and Pethia nigrofasciata Mitochondrial Genomes.

Smiliogastrinae are recognized for their high nutritional and ornamental value. In this study, we employed high-throughput sequencing technology to acquire the complete mitochondrial genome sequences of Dawkinsia filamentosa and Pethia nigrofasciata. The gene composition and arrangement order in these species were similar to those of typical vertebrates, comprising 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, and 1 non-coding region. The mitochondrial genomes of D. filamentosa and P. nigrofasciata measure 16,598 and 16,948 bp, respectively. Both D. filamentosa and P. nigrofasciata exhibit a significant preference for AT bases and an anti-G bias. Notably, the AT and GC skew values of the ND6 gene fluctuated markedly, suggesting that the selection and mutation pressures on this gene may differ from those affecting other genes. Phylogenetic analysis, based on the complete mitochondrial genomes of 23 Cyprinidae fishes, revealed that D. filamentosa is closely related to the sister group comprising Dawkinsia denisonii and Sahyadria chalakkudiensis. Similarly, P. nigrofasciata forms a sister group with Pethia ticto and Pethia stoliczkana.

The tri-flow adaptiveness of codes in major evolutionary transitions.

Life codes increase in both number and variety with biological complexity. Although our knowledge of codes is constantly expanding, the evolutionary progression of organic, neural, and cultural codes in response to selection pressure remains poorly understood. Greater clarification of the selective mechanisms is achieved by investigating how major evolutionary transitions reduce spatiotemporal and energetic constraints on transmitting heritable code to offspring. Evolution toward less constrained flows is integral to enduring flow architecture everywhere, in both engineered and natural flow systems. Beginning approximately 4 billion years ago, the most basic level for transmitting genetic material to offspring was initiated by protocell division. Evidence from ribosomes suggests that protocells transmitted comma-free or circular codes, preceding the evolution of standard genetic code. This rudimentary information flow within protocells is likely to have first emerged within the geo-energetic and geospatial constraints of hydrothermal vents. A broad-gauged hypothesis is that major evolutionary transitions overcame such constraints with tri-flow adaptations. The interconnected triple flows incorporated energy-converting, spatiotemporal, and code-based informational dynamics. Such tri-flow adaptations stacked sequence splicing code on top of protein-DNA recognition code in eukaryotes, prefiguring the transition to sexual reproduction. Sex overcame the spatiotemporal-energetic constraints of binary fission with further code stacking. Examples are tubulin code and transcription initiation code in vertebrates. In a later evolutionary transition, language reduced metabolic-spatiotemporal constraints on inheritance by stacking phonetic, phonological, and orthographic codes. In organisms that reproduce sexually, each major evolutionary transition is shown to be a tri-flow adaptation that adds new levels of code-based informational exchange. Evolving biological complexity is also shown to increase the nongenetic transmissibility of code.

Anciently duplicated genes continuously recruited to heart expression in vertebrate evolution are associated with heart chamber increase.

Although gene/genome duplications in the early stage of vertebrates have been thought to provide major resources of raw genetic materials for evolutionary innovations, it is unclear whether they continuously contribute to the evolution of morphological complexity during the course of vertebrate evolution, such as the evolution from two heart chambers (fishes) to four heart chambers (mammals and birds). We addressed this issue by our heart RNA-Seq experiments combined with published data, using 13 vertebrates and one invertebrate (sea squirt, as an outgroup). Our evolutionary transcriptome analysis showed that number of ancient paralogous genes expressed in heart tends to increase with the increase of heart chamber number along the vertebrate phylogeny, in spite that most of them were duplicated at the time near to the origin of vertebrates or even more ancient. Moreover, those paralogs expressed in heart exert considerably different functions from heart-expressed singletons: the former are functionally enriched in cardiac muscle and muscle contraction-related categories, whereas the latter play more basic functions of energy generation like aerobic respiration. These findings together support the notion that recruiting anciently paralogous genes that are expressed in heart is associated with the increase of chamber number in vertebrate evolution.

A CTCF-dependent mechanism underlies the Hox timer: relation to a segmented body plan.

During gastrulation, Hox genes are activated in a time-sequence that follows the order of the genes along their clusters. This property, which is observed in all animals that develop following a progressive rostral-to-caudal morphogenesis, is associated with changes in the chromatin structure and epigenetic profiles of Hox clusters, suggesting a process at least partly based on sequential gene accessibility. Here, we discuss recent work on this issue, as well as a possible mechanism based on the surprising conservation in both the distribution and orientation of CTCF sites inside vertebrate Hox clusters.

Development and validation of a universal primer pair for the taxonomic and phylogenetic studies of vertebrates.

BACKGROUND: Recent studies in the field of molecular identification have described 16S rRNA gene as a highly informative fragment of mitochondrial DNA for species discrimination. This study presents a newly developed universal primer pair yielding an approximately 350 bp fragment of mitochondrial 16S rRNA, variable enough to encompass and identify all vertebrate classes. METHODS AND RESULTS: The primers were designed by aligning and analyzing over 1500 16S rRNA sequences downloaded from the NCBI nucleotide database. A total of 93 vertebrate species, spanning 27 orders and 55 families, were PCR-amplified to validate the primers. All the target species were successfully amplified and identified when aligned with reference sequences from the NCBI nucleotide database. Using the Kimura 2-parameter model, low intra-species genetic divergence of the target region was observed - from 0 to 4.63%, whereas relatively higher inter-species genetic divergence was observed, ranging from 4.88% to 69.81%. Moreover, the newly developed primers were successfully applied to a direct PCR protocol, making the workflow very cost-effective, time-saving and less laborious in comparison to conventional PCR. CONCLUSIONS: The short length, high variability and conserved priming sites of the target fragment across all vertebrate species make it a highly desirable marker for species identification and discrimination.

Sea lamprey enlightens the origin of the coupling of retinoic acid signaling to vertebrate hindbrain segmentation.

Retinoic acid (RA) is involved in antero-posterior patterning of the chordate body axis and, in jawed vertebrates, has been shown to play a major role at multiple levels of the gene regulatory network (GRN) regulating hindbrain segmentation. Knowing when and how RA became coupled to the core hindbrain GRN is important for understanding how ancient signaling pathways and patterning genes can evolve and generate diversity. Hence, we investigated the link between RA signaling and hindbrain segmentation in the sea lamprey Petromyzon marinus, an important jawless vertebrate model providing clues to decipher ancestral vertebrate features. Combining genomics, gene expression, and functional analyses of major components involved in RA synthesis (Aldh1as) and degradation (Cyp26s), we demonstrate that RA signaling is coupled to hindbrain segmentation in lamprey. Thus, the link between RA signaling and hindbrain segmentation is a pan vertebrate feature of the hindbrain and likely evolved at the base of vertebrates.

Global conservation status of the jawed vertebrate Tree of Life.

Human-driven extinction threatens entire lineages across the Tree of Life. Here we assess the conservation status of jawed vertebrate evolutionary history, using three policy-relevant approaches. First, we calculate an index of threat to overall evolutionary history, showing that we expect to lose 86-150 billion years (11-19%) of jawed vertebrate evolutionary history over the next 50-500 years. Second, we rank jawed vertebrate species by their EDGE scores to identify the highest priorities for species-focused conservation of evolutionary history, finding that chondrichthyans, ray-finned fish and testudines rank highest of all jawed vertebrates. Third, we assess the conservation status of jawed vertebrate families. We found that species within monotypic families are more likely to be threatened and more likely to be in decline than other species. We provide a baseline for the status of families at risk of extinction to catalyse conservation action. This work continues a trend of highlighting neglected groups-such as testudines, crocodylians, amphibians and chondrichthyans-as conservation priorities from a phylogenetic perspective.

Unstructured linker regions play a role in the differential splicing activities of paralogous RNA binding proteins PTBP1 and PTBP2.

RNA Binding Proteins regulate, in part, alternative pre-mRNA splicing and, in turn, gene expression patterns. Polypyrimidine tract binding proteins PTBP1 and PTBP2 are paralogous RNA binding proteins sharing 74% amino acid sequence identity. Both proteins contain four structured RNA-recognition motifs (RRMs) connected by linker regions and an N-terminal region. Despite their similarities, the paralogs have distinct tissue-specific expression patterns and can regulate discrete sets of target exons. How two highly structurally similar proteins can exert different splicing outcomes is not well understood. Previous studies revealed that PTBP2 is post-translationally phosphorylated in the unstructured N-terminal, Linker 1, and Linker 2 regions that share less sequence identity with PTBP1 signifying a role for these regions in dictating the paralog's distinct splicing activities. To this end, we conducted bioinformatics analysis to determine the evolutionary conservation of RRMs versus linker regions in PTBP1 and PTBP2 across species. To determine the role of PTBP2 unstructured regions in splicing activity, we created hybrid PTBP1-PTBP2 constructs that had counterpart PTBP1 regions swapped to an otherwise PTBP2 protein and assayed on differentially regulated exons. We also conducted molecular dynamics studies to investigate how negative charges introduced by phosphorylation in PTBP2 unstructured regions can alter their physical properties. Collectively, results from our studies reveal an important role for PTBP2 unstructured regions and suggest a role for phosphorylation in the differential splicing activities of the paralogs on certain regulated exons.

Diversity and evolution of the vertebrate chemoreceptor gene repertoire.

Chemoreception - the ability to smell and taste - is an essential sensory modality of most animals. The number and type of chemical stimuli that animals can perceive depends primarily on the diversity of chemoreceptors they possess and express. In vertebrates, six families of G protein-coupled receptors form the core of their chemosensory system, the olfactory/pheromone receptor gene families OR, TAAR, V1R and V2R, and the taste receptors T1R and T2R. Here, we study the vertebrate chemoreceptor gene repertoire and its evolutionary history. Through the examination of 1,527 vertebrate genomes, we uncover substantial differences in the number and composition of chemoreceptors across vertebrates. We show that the chemoreceptor gene families are co-evolving, highly dynamic, and characterized by lineage-specific expansions (for example, OR in tetrapods; TAAR, T1R in teleosts; V1R in mammals; V2R, T2R in amphibians) and losses. Overall, amphibians, followed by mammals, are the vertebrate clades with the largest chemoreceptor repertoires. While marine tetrapods feature a convergent reduction of chemoreceptor numbers, the number of OR genes correlates with habitat in mammals and birds and with migratory behavior in birds, and the taste receptor repertoire correlates with diet in mammals and with aquatic environment in fish.

Reference transcriptome assembly of a protogynous sex change fish, harlequin sandsmelt (Parapercis pulchella).

The harlequin sandsmelt (Parapercis pulchella) is a female-to-male sex change fish in which functional females possess ovotestes that consist of both ovarian and testicular tissues. These features indicate that this species could be an excellent model for studying the flexibility of sex differentiation in vertebrates. However, genetic resources in this species have so far been limited. Therefore, in this study, the reference transcriptome of this fish was constructed through RNA-sequencing, de novo transcriptome assembly, superTranscripts construction, and functional annotations. To obtain as many genes as possible, RNA was extracted from various tissues (brains, gills, hearts, livers, guts, and gonads) and various sexual stages (females, individuals during sex change, and males) and then subjected to sequencing and downstream analyses. As a result, 91,884 representative transcripts with 32,627 protein-coding sequences were generated. 72.2% of protein-coding sequences (23,566 sequences) were functionally annotated. Also, our analysis shows that the superTranscripts method effectively removes redundant sequences from raw-assembled data compared with other strategies. The resultant dataset is a valuable resource for future molecular developmental studies on sex change in P. pulchella.